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ABSTRACT 



THE KINEMATICS AND PHYSICAL CONDITIONS OF THE IONIZED 
GAS IN NGC 4593 FROM CHANDRA HIGH-ENERGY GRATING 

SPECTROSCOPY 

ff^ ■ Barry McKernan^, Tahir Yaqoob^'^, Ian M. George^'^, T. J. Turner'^'^ 

o ■ 
o : 

^ ■ 

. We observed the Seyfert 1 galaxy NGC 4593 with the Chandra high energy trans- 

! mission gratings and present a detailed analysis of the soft X-ray spectrum. We measure 

I strong absorption lines from He-like O, Ne, Mg, Si, H-like N, O, Ne, Mg, Si and highly 

I ionized Fe xix-xxv. The weighted mean of the offset velocity of the strongest absorp- 

^T) ' tion profiles is —140 it 35 km s~^. However the individual profiles are consistent with 

■ the systemic velocity of NGC 4593 and many profiles hint at the presence of either 
t:;^ ■ multiple kinematic components or blending. Only the N vii Lya (A24.781) , O viii Lya 
O , (A18.969) and Mg xii Lya (A8.421) lines appear to be marginally resolved. We identify 
Q . a spectral feature at ~ 0.707 keV with a neutral Fe L edge, which might suggest that 
^ I there is dust along the line-of-sight to NGC 4593 although, this is not the only inter- 

C^l pretation of this feature. A search for neutral O absorption which would reasonably be 

O ' expected from dust absorption is complicated by contamination of the Chandra ACIS 

' CCDs. Neutral Si absorption, which might also be expected from absorption due to 

■ dust is present (though not significantly) in the form of a weak neutral Si edge. The 
J> . neutral Si column (< 4 x 10^'^ cm~^) corresponding to the Si edge is consistent with 

^ , the neutral Fe column (~ 1.5 x 10^^ cm^^) from the Fe L edge. We also detect, at 

^ ! marginal signficance, N vii Lya (A24.781) and O vii (r) Is^ — ls2p (A21.602) absorp- 

tion at z ~ 0, due to a hot medium in our Local Group. The soft X-ray spectrum of 
NGC 4593 is adequately described by a simple, single-zone photoionized absorber with 
an equivalent Hydrogen column density of 5.37~^qyq x 10^^ cm~^ and an ionization pa- 
rameter of log^ = 2.52+°;°^ ergs cm s although there remain some features which are 
not identified. Although the photoionized gas almost certainly is comprised of matter in 
more than one ionization state and may consist of several kinematic components, data 
with better signal-to-noise ratio and better spectral resolution are required to justify a 
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more complex model. Finally, in emission we detect only weak forbidden lines of [Ne ix] 
and [O Vll] . 

Subject headings: galaxies: active - galaxies: individual (NGC 4593 ) - galaxies: Seyfert 
- techniques: spectroscopic - ultraviolet: galaxies - X-rays: galaxies - X-rays: galax- 
iesgalaxies: active - galaxies: individual (NGC 4593 ) - galaxies: Seyfert - techniques: 
spectroscopic - X-rays: galaxies 

1. Introduction 

The launch of Chandra and XMM-Newton began a new era in the study of X-ray photoionized 
circumnuclear gas. High resolution X-ray spectroscopy with Chandra now allows the gas kine- 
matics to be studied seriously for the first time, and the detection of individual absorption and 
emission lines can now place very strong constraints on the ionization structure of the warm, or 
partially ionized, X-ray absorbing gas in type 1 active galactic nuclei (AGNs)(e.g. (Lee et al. 2001); 
(Sako et al. 2001a); (Kaastra et al. 2002); (Kaspi et al. 2002); (Yaqoob et al. 2003)). 

Here we present the results of a Chandra High-Energy Transmission Grating (or HETGS - 
Markert, et al. 1995) observation of NGC 4593 , which was simultaneous with an i^XTE' observation. 
NGC 4593 {z = 0.0083^, De Vaucouleurs et al. 1991) is a type I AGN that lies within a barred- 
spiral host galaxy. NGC 4593 is also relatively luminous (L2-10 keV typically 8.6 x 10^^ ergs s~^, 
converted from the luminosity given in Reynolds (1997) using Hq = 70 km s~^ Mpc~^ and qo = 
which we use throughout this paper). In the present paper we focus on the soft X-ray spectroscopy 
results. 

The paper is organized as follows. In §2 we present the data and describe the analysis tech- 
niques. In §3 we discuss gross features of the X-ray spectrum, including the intrinsic continuum 
form and interpret the data in the context of historical, lower spectral resolution CCD data. In 
§4 we qualitatively discuss the discrete X-ray spectral features before describing detailed spectral 
modeling. In §5 we describe in detail the modeling of the X-ray spectrum using the photoionization 
code, XSTAR. In §6 we compare our results with those from other X-ray observations of NGC 4593 
. In §7 we discuss the possibility of a dusty warm absorber in NGC 4593 and we investigate the 
relativistically broadened disk-line scenario (Branduardi-Raymont et al. 2001; Sako et al. 2001b; 
Mason et al. 2002) as an alternative to the dusty warm absorber model. Finally, in §8 we summarize 
our results and in §9 briefly state our conclusions. 



^Deduced from observations of the 21cm H I line. Uncertainty in this measurement (±6 km s~^) is negligible 
compared to the systematic uncertainty in the Chandra MEG energy scale (±67 km s~^ at 0.5 keV) and also much 
less than the best MEG spectral resolution (~ 280 km s~^ at 0.5 keV). Optical observations suggest a range of 
z = 0.0080 — 0.0090, corresponding to a range of —90 km s~^ to +180 km s~^. However, optical observations may 
be confused by AGN outflow and so are less accurate indicators of source redshift. 
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2. Observations and Data 

We observed NGC 4593 with Chandra (simultaneously with RXTE ) on 2001 June 29 for a 
duration of ~ 79 ks, beginning at UT 00:54:49. The Chandra data were reprocessed using ciao 
2.1.3 and CALDB version 2.7, according to recipes described in ciao 2.1.3 threads^. The RXTE 

PC A data were reduced using methods described in Weaver, Krolik, & Pier (1998), except improved 
background subtraction was utilized (released 2002, February) and a later version of the spectral 
response matrix (v 8.0) was used®. Only data from layer 1 and from PCUs 1, 2, 3 and 4 were used 
(PCUO has had technical problems in the later part of the mission). 

For Chandra the instrument used in the focal plane of the High Resolution Mirror Assembly was 
the HETGS . The HETGS consists of two grating assemblies, a High-Energy Grating (HEG) and 
a Medium-Energy Grating (MEG). Only the summed, negative and positive, first-order Chandra 
grating spectra were used in our analyis. For further discussion of the HETGS and our analysis 
procedure here, see Yaqoob et al. (2003). The mean Chandra total HEG and MEG count rates 
were 0.1965 ± 0.0017 and 0.4743 ± 0.0021 cts/s respectively. HEG, MEG and PGA spectra were 
extracted over the entire on-time for each instrument. This resulted in net exposure times of 78,889 
s for HEG and MEG, and between ~ 14 ks and ~ 70 ks for each PGU. The background subtracted 
spectra from the different PCUs were co-added for analysis and the resulting mean count rate was 
4.8 ct/s over the entire PCA band. The excess variances on the hardness ratios 1.0-2.5 keV/(0.6- 
1.0 keV) and 2.5-7.0 keV/(0.6-1.0 keV) for Chandra data binned at 0.161 were 0.0011 ± 0.0012 
and 0.0014 it 0.0007 respectively. Thus, the data are consistent with no spectral variability in soft 
X-rays during the observation, but there may be some small spectral variation in hard X-rays. The 
RXTE observation extended ~ 20 ks beyond the Chandra observation but the mean count rate in 
that portion of the observation was compatible with the mean of the strictly overlapping portion. 
Therefore, for the sake of improved statistics, the PCA spectrum included ~ 10 ks exposure beyond 
the Chandra data. 

For Chandra HETGS we made effective area files (ARF, or ancillary response file), photon 
spectra and counts spectra following the method of Yaqoob et al. (2003). Since the MEG soft 
X-ray response is much better than the HEG (whose bandpass only extends down to ~ 0.8 keV) 
we will use the MEG as the primary instrument but refer to the HEG for confirmation of features 
in the overlapping bandpass and for constraining the continuum. The MEG spectral resolution 
corresponds to FWHM velocities of ~ 280,560, and 3,560 km s~^ at observed energies of 0.5, 1.0, 
and 6.4 keV respectively. 

We did not subtract detector or X-ray background since it is such a small fraction of the 
observed counts (< 0.5% for the MEG). We treated the statistical errors on both the photon and 
counts spectra with particular care since the lowest and highest energies of interest can be in the 



^http: / / asc.harvard.edu/ciao2. 1 /documents-threads. html 

^Details of both in http://heasarc.gsfc.nasa.gov/docs/xte/xhp_proc_analysis.html 
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Poisson regime, with spectral bins often containing a few, or even zero counts. We assign statistical 
upper and lower errors of 1.0+ ^/(iV + 0.75) and iV(1.0- [1.0- (l/9iV) - (l/3y(iV))]^) respectively 
(from Gehrels, 1986) on the number of photons, N, in a given spectral bin. Statistical errors were 
always computed on numbers of photons in final bins, according to the above prescription for 
plotting purposes. When fitting the Chandra data, we used the C-statistic (Cash, 1976) for finding 
the best-fitting model parameters, and quote 90% confidence, one-parameter statistical errors unless 
otherwise stated. The C— statistic minimization algorithm is inherently Poissonian and so makes 
no use of the errors on the counts in the spectral bins described above. 

For the HEG data, the signal-to-noise ratio per 0.02A bin is > 1 blueward of ~ 15A and as 
high as ~ 10 in the 6 — 7 A range. For the MEG data, in the ^ 2 — 25A band which we examine 
here, the signal-to-noise ratio per 0.02A bin ranged from ~ 1 at ~ 22.5A to a maximum of ~ 14 
at ~ lA. The systematic uncertainty in the energy scale is currently believed to be 0.0028^ and 
0.00551 for the HEG and MEG respectively^. For the MEG, at 0.5 keV, 1 keV and 6.4 keV this 
corresponds to velocity offsets of ~ 67, 133, and 852 km s~^ respectively, and about half of these 
values for the HEG. 

In addition, there has been a continuous degradation of the quantum efficiency of Chandra 
ACIS due to molecular contamination. The degradation includes an excess absorption due to 
neutral O ^. When fitting photoionization models to the data we will investigate the effect of this 
degradation on model parameters. There is no satisfactory correction available for the gratings but 
we will estimate worst case effects using a pure ACIS correction, since a grating in front of the 
CCDs will not make the CCD degradation worse. Indeed a grating may improve the situation by 
adding extra energy information. We will show in §5.3 that omitting the correction only affects the 
inferred intrinsic continuum and does not affect the important physical parameters of absorption 
models. 



3. Overall Spectrum and Preliminary Spectral Fitting 

We used XSPEC vll.2.0 for spectral fitting to the HEG and MEG spectra in the 0.8-5 keV 
and 0.5-5 keV bands respectively. These energy bands will be used in all the spectral fitting in 
the present paper, in which we concentrate on features in the soft X-ray spectrum (less than ~ 2 
keV). The harder spectrum out to ~ 15 keV, including the Fe-K emission line and Compton- 
refiection continuum (using simultaneous RXTE data in addition to the Chandra data) will be 
discussed elsewhere. However, we will utilize a portion of the RXTE PCA spectrum to constrain 
the continuum for fitting the soft X-ray Chandra data and use the statistic for minimization. 

Since we convolve all models through the instrument response before comparing with the data. 



'^http://space.mit.edu/CXC/calib/hetgcal.html 

*http://cxc.harvard.edu/cal/Links/Acis/acis/CaLprods/qeDeg/index.html 
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all our model line-widths axe intrinsic values and do not need to be corrected for the instrument 
broadening. This contrasts with the analysis method in some of the literature on results from 
grating data in which models are fitted directly to the data and then corrected for instrument 
broadening (e.g. Kaspi et al. 2002). First, we show how the 0.5-5 keV MEG data compare to 
a simple model consisting only of a single power law and Galactic absorption. For the latter, we 
use a value of 1.97 x lO^'^cm"^ (Elvis et al. 1989) throughout this work. Later in this paper we 
show that the ACTS degradation discussed above does not effect our results. Therefore we do not 
let the Galactic column vary in spectral fitting since its systematic error is much less than the 
effects of ACIS degradation. Fig. 1 (a) shows the MEG photon spectrum of NGC 4593 binned 
at 0.32^ (which is comparable with ASCA resolution) with the model overlaid. The power-law 
index was allowed to vary in this fit. Fig. 1 (b) shows the ratio of the data to this model. This 
simple power law is clearly a poor fit to the data, with considerable complexity apparent, both in 
the continuum, and in terms of discrete absorption features. Fig. 1 (b) clearly shows a soft X-ray 
excess, rising up from the hard power law (modified by Galactic absorption), below 0.7 keV. We 
caution that there may still be broad residual uncertainties of ~ 30% or more in the calibration of 
the MEG effective area at the lowest energies ^. The instrumental effective area is very small at 
these energies, although residuals at < 0.6 keV probably coincide with the complex instrumental 
O I edge which is spread over ~ 10 eV around ~ 0.53 keV. 

Observations of NGC 4593 in 1994 with ASCA , as summarized in Reynolds (1997) and George 
et al. (1998) also indicated a soft excess present below ~ 0.7 keV, relative to an extrapolation of 

the best-fitting hard X-ray power-law. A soft excess has also been observed by XMM-Newton and 
Chandra LETGS (Kaastra et al. 2001), as well as by BeppoSAX (Guainazzi et al. 1999). This 
would suggest that at least part of the soft excess we observed with Chandra HETGS is real. The 
low-energy degradation of the CCDs mentioned in §2 is in fact in the sense that any soft excess 
will actually be larger when the correction is taken into account. 

The origin of the soft excess is unclear. It may be the tail end of some kind of soft thermal 
emission (or Comptonized soft thermal emission), possibly from an accretion disk (e.g. see Piro, 
Matt, k, Ricci (1997), and references therein). Any relativistically-broadened O Vlll Lya (A18.969) 
emission that is present could also contribute to the soft excess (e.g. Branduardi-Raymont et al. 
2001; Turner et al. 2001). Since the soft excess appears only in the 0.5-1 keV band of our data, 
we do not have enough information to constrain its origin, and sophisticated modeling of it is not 
warranted. Therefore, in the remainder of this paper we model the 0.5-5 keV intrinsic continuum 
with a broken power law. When the data are modeled with this continuum, there is still considerable 
structure in the soft X-ray spectrum, characteristic of the presence of a photoionized, or 'warm' 
absorber. 

Now, the warm absorber in NGC 4593 and other AGN when observed with CCDs (such as 
those on ASCA ) has often been modeled simply with absorption edges due to O Vll and O Vlll 
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since the lower CCD spectral resolution did not usually warrant more sophisticated models. In 
order to directly compare the new Chandra data with previous observations we modeled the MEG 
data for NGC 4593 with two absorption edges, Galactic absorption, and the intrinsic continuum 
described above. Here we used the RXTE PCA data to help constrain the continuum. First, we 
modeled the HEG, MEG, and PCA data simultaneously, the latter data only covering the 3-8 
keV band in order to avoid the Compton-reflection continuum. We had to include an additional 
Gaussian model component (with center energy, width and intensity allowed to vary) to model 
the Fe-K line. Since we are not concerned with the Fc-K line in this paper we then take the hard 
power-law index (F = 1.794^QQj^g) from this fit and use T = 1.794 in all subsequent model-fitting 
in this paper. If we allow F to vary in our best model (two absorption edges, a warm absorber 
and a broken power- law , see §8 below) fit to the Chandra MEG data binned at 0.16A, we find a 
best-fitting value of F = 1.803 it 0.017 which is consistent with the fixed value above, so fixing F 
does not affect our results or conclusions. Since we will only be using Chandra data below 5 kcV 
from this point on, we do not have to worry about the spectral fits becoming unstable due to too 
much interplay between the two parts of the broken power-law index and the Fe-K line model. This 
will be dealt with in future work. 

In order to measure the absorption-edge parameters, we then modeled the 0.5-5 kcV MEG 
data only, fixing the hard photon index at 1.794, and allowed the energies and optical depths (at 
threshold) of the two absorption edges to vary. The best-fitting model overlaid on the MEG photon 
spectrum is shown in Fig. 2 . Although a spectrum with bin size 0.04A was used for the spectral 
fitting, the plot shows the data binned at 0.32A since this is closer to the resolution of CCD data. 
The best-fitting parameters obtained from this model were 0.708 it 0.003 keV and 0.868^0050 
for the edge threshold energies. For the the optical depths at threshold we obtained r = 0.18 it 0.04 
and 0.1410:03 respectively. All quantities refer to the source rest frame. The best-fitting soft X-ray 
photon index was 2.26 ± 0.05 and the break energy was 1.031q q3 

The edge at 0.868 keV is in good agreement with the expected energy of the O viii edge at 0.871 
keV. However, the edge at 0.708 keV is redshifted by ~ 12, 600 km s~^ from the expected energy of 
the O VII edge (0.739 keV). Rather, this edge energy agrees extremely well with the expected energy 
of a neutral Fc L3-cdge at 0.707 keV (as observed, for example, by Lee et al. (2001) in the warm 
absorber of MCG —6—30—15 ). Furthermore, adding a third edge to the continuum model, with a 
threshold energy fixed at that expected for O Vll , does not improve the fit statistic significantly. 
The threshold optical depth of the O vii edge in this case is < 0.03 at the one-parameter, 90% 
confidence level. We note that since the spectrum is so complex around the regions of the O Vll 
and O VIII edges, optical depths derived from simple models should be interpreted with caution. 
We can compare the threshold optical depths with those measured from the 1994 January 4 ASCA 
observation of NGC 4593 (Reynolds 1997): 0.26 ± 0.04 (ascribed to O vii ) and 0.09l:[|;[J^ (O Vlll 
). The Chandra and ASCA O Vlll edge depths arc in marginal agreement but the O Vll depths 
are completely incompatible. It seems that, due to the limited energy resolution of CCD data, 
the Fe-L3 edge was misidentified with an O Vll edge, not only in NGC 4593 but possibly also for 



- 7- 



other AGN. The high resolution of Chandra HETGS indicates that this edge cannot be due to 
O VII , and previous ASCA data from AGN needs to be re-examined in this context. The optical 
depth that we measure for the Fe-L3 edge is in marginal agreement with that ascribed to the 
O VII edge by Reynolds (1997). The ASCA SIS (4-CCD mode) observation of NGC 4593 had an 
energy resolution (FWHM) of ~ 140 eV at ~ 1.0 keV at the time of the observation, and therefore 
could not distinguish between an O Vll edge and a neutral Fe edge. However, an Fe-L edge for 
the feature at 0.7 keV is not the only interpretation (see §7). High resolution observations with 
Chandra are also hinting at X-ray absorption fine structure (XAF) around absorption edges (Lee 
et al. 2002). The apparent Fe-L edge in our data extends over ~ 5 eV, which is ~ 3 times the 
FWHM resolution of Chandra at this energy. The S/N of our data and the FWHM resolution of 
Chandra are not adequate to measure the XAF structure around the Fe-L edge, which in any case, 
will be complicated by O Vll Is^ — Isnp (n > 4) absorption features. 



4. Spectral Features and Kinematics 

4.1. Absorption and Emission Lines 

In order to best illustrate different characteristics of the spectrum wc display the spectral 
data in several different ways. In Fig. 3 we show the MEG photon spectrum between 0.47-2.1 
keV, as a function of observed energy. The data and the models in Fig. 3 are all absorbed by the 
Galactic column density. An important point to realize from Fig. 3 is that in some regions the 
spectrum is so complex that it is difficult to distinguish absorption from emission. For example, 
two emission features may be present between 1.0-1.5 keV according to the data/model ratio plot 
of Fig. 2 . However, from Fig. 3 it is clear that the peak flux of both features coincides with the 
inferred intrinsic continuum, implying that the apparent feature is simply due to the presence of 
broad absorption complexes on either side of it. In addition there is a sharp variation in the MEG 
effective area (at ~ 0.85 keV) and any uncertainties in the response could be manifested as an 
apparent spectral feature (see §5 for further discussion). 

In Fig. 4 and Fig. 5 we show the MEG spectrum, this time as a function of wavelength in the 
source rest-frame, redward of 2A. These spectra have bin sizes of 0.02A, approximately the FWHM 
MEG spectral resolution (0.023A). In Fig. 4 we overlay the Lyman series (blue) for H-like N, O, 
Ne, Mg, Si, S, Ar and the corresponding He-like triplets (red) that lie in the 2-25A range. In Fig. 5 
we overlay the He-like resonance series lines (blue) to the n = 1 level for the same elements. In 
addition, in Fig. 5 we indicate L-shell transitions of Fe XIX (red) in the 2 — 25A range^° to illustrate 
just some of the blending of spectral features that can arise due to absorption by highly ionized Fe. 

Referring to Figs. 3 to 5, we see that the strongest, clear discrete spectral features detected 
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are absorption lines due to O Vlll Lya (A18.969) , O Vlll Ly/? (A16.006) , Ne IX (r) Is^ - ls2p 
(A13.447), Ne ix (r) Is^ - ls3p (All. 549), Ne X Lya (A12.134) , Mg Xll Lya (A8.421) , Si xm 
(r) Is^ — ls2p (A6.648). Also included amongst the strongest discrete features are transitions (and 
blends of transitions) which we identify with highly ionized Fe, typically L-shell transitions in Fe xx 
-Fe xxv . Weaker absorption lines are also apparent from the spectra such as N Vli Lya (A24.781) 
, O VII (r) Is^ — ls2p (A21.602) and higher-order Lyman and (He-like) resonance absorption lines 
of O and Ne (there is evidence for transitions up to n = 4 in each case; see Figs. 4 and 5). Some 
of these weaker features may be blended with Fe XXII - Fe XXIII transitions. Higher signal-to-noise 
X-ray spectroscopy (e.g. of NGC 3783, Kaspi et al. 2002) has shown that blending of absorption 
features considerably complicates the spectrum. In particular, transitions of Fe xxiii could be 
blended with the absorption troughs at ~ 10.65 and 11. OA. Similarly, transitions of Fe xx and 
Fe XXII could be blended with troughs at 9.15A, 12. 8A, and 11.75A respectively. We note that such 
a wide range of observed ionization states (from N VII to Fe xxv ) suggests that there is more than 
one ionization component of the warm absorber. We shall return to this point in the discussion of 
photoionization modeling in §5. 

We also identify absorption features due to O Vll (r) Is''^ — ls2p (A21.602) and N VII Lya 
(A24.781) in the observed frame, at z = 0, indicating absorption along the line of sight due to a 
hot medium in our Galaxy or local group (see also Figs. 3, 4 &; 5). The separations of the lines at 
2; = and z = 0.0083 is 2490 km s~^, which is easily resolved by Chandra . Strong local absorption 
by O VII (r) Is^ - ls2p (A21.602) and Ne IX (r) Is^ - ls2p (A13.447) and weaker local absorption 
by O VIII Lya (A18.969) and O Vli (r) Is^ - ls3p (A18.628) has also b een observed by Nicastro et 
al. (2002) in the direction of the blazar PKS 2155-304. Nicastro et al. (2002) suggest that the local 
absorption is produced in a low density intergalactic plasma, collapsing towards our Galaxy. It is 
likely that the 2; = absorption features O vii (r) Is^ - ls2p (A21.602) and N vii Lya (A24.781) in 
NGC 4593 are due to the same gas. The strengths of the lines, although only marginally significant, 
are comparable with the EW~ 0.2 — 0.3 eV detected by Nicastro et al. (2002). 

Apart from Fe Ka (which will be discussed elsewhere) , we also identify weak [O Vll] and [Ne ix] 
emission features. Table 1 shows that [O Vll] and [Ne ix] are redshifted but the errors are consistent 
with systemic velocity. The corresponding intercombination and resonance lines of the O Vll and 
Ne IX He-like triplets are not detected. The relative prominence of forbidden lines compared to the 
intercombination and resonance lines indicates an origin in a photoionized plasma (e.g. Porquet & 
Dubau 2000). Fig. 6 shows close-ups of the wavelength regions spanning the O vii and Ne ix He- like 
triplets (with the positions of the forbidden (f ) , intercombination (i) and resonance (r) lines marked 
in red). The ratio of the forbidden and intercombination line fluxes (R=f/i) constrains the electron 
density (e.g. see Porquet & Dubau 2000). An estimate of R from the flux at the i,f wavelengths 
in Fig. 6 yields R<3 for Ne ix and R<10 for O Vll . These ratios imply rig < 2 x 10^^ cm~^ and 
ng < 8 X 10"'^'^ cm^'"^ respectively, for a relative ionic abundance of the H-likc and He-like ions of 
< 10'^ (see Fig. 9 of Porquet & Dubau 2000). If blending is significant in either of these triplets, then 
the measured line intensities would be upper limits only, thereby yielding a potentially significant 
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systematic error to our estimates of R and therefore Ue- 

The presence of broad, bound-free absorption features and many narrow features which may be 
blended, makes it impossible to define an observed continuum for much of the spectrum, and there- 
fore it is difficult to measure meaningful equivalent widths of even the strongest narrow absorption 
features. Rather than try to measure equivalent widths of all candidate features, and compare 
with model predictions, we will take the physically more direct approach of fitting photoionization 
models to the data (§5). First we give the measured energies and equivalent widths of the strongest 
absorption lines. The measurements were made using the broken power law plus two-edge model as 
the baseline spectral model (see §3) and Gaussians to model the absorption lines. The continuum 
and edge parameters were frozen at the best-fitting values given in §3, and the data were fitted 
over a few hundred eV centered on the feature of interest. Allowing the edge parameters to vary 
did not change our results. The fitted range of the data was extended down to 0.47 keV in order 
to measure the parameters of N Vll Lya (A24.781) . Initially the intrinsic width of each Gaussian 
was frozen at less than the instrumental resolution but then allowed to vary in order to investigate 
whether the lines were resolved. The results for the measured energies, apparent velocity shifts 
relative to systemic, equivalent widths, intrinsic widths and the improvement in the fit-statistic, of 
these strongest absorption features are given in Table 1. Note that the quoted equivalent widths 
are relative to the intrinsic continuum (corrected for Galactic absorption). Again, we caution that 
the measured equivalents widths are subject to uncertainties in the continuum, line-blending, and 
any complexity in the line profile. It is instructive to compare line strengths with a physical model, 
so also in Table 1 are the predicted equivalent widths from the best-fitting photoionization model 
discussed in S5. 



4.2. Absorption-Line Velocity Profiles 

Table 1 shows that most of the measured absorption lines are blucshifted, with Gaussian- fitted 
centroids ranging from near systemic velocity to ~ —310 km s^"*^. The exceptions are Ne IX (r) 
Is^ — ls2p (A13.447) and Fe xx (A12.817). The former is redshifted, but the measurement is 
consistent with systemic velocity and the latter is likely blended with nearby Fe XX transitions 
(see below). Some of the velocity profiles of the absorption features in Table 1 appear to be more 
complex than Gaussian, as can be seen from an inspection of the profiles in Fig. 7 and Fig. 8. These 
velocity profiles were obtained from combined HEG and MEG data except for those transitions at 
> 17 A, which were determined using MEG data only, since they are out of the HEG bandpass. The 
profiles are centered with zero at the systemic velocity of NGC 4593 , with positive and negative 
velocities corresponding to redshifts and blueshifts respectively. Superimposed on the data are 
velocity spectra from the best-fitting photoionization model which will be discussed in detail in §5. 

Fig. 7 (a) shows the velocity profile centered on the O Vlll Lya (A18.969) absorption transition 
in the NGC 4593 rest-frame. A broad absorption trough can be seen between ~ —600 km s"^ and 
systemic, centered around ~ —300 km s""^. Fig. 7 (b) is centered around the Ne X Lya (A12.134) 
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transition and also exhibits a broad absorption profile from ~ —300 km s~^ to ^ +100 km s~^. 
Fig. 7 (c) shows the velocity profile for Mg XII Lya (A8.421) , which exhibits a broad absorption 
profile from ^ —500 km s^^ to +300 km s^-*^. Fig. 7 (d) shows the velocity profile for Si xiv Lya 
(A6.182) . The Si xiv Lya (A6.182) profile appears to be broader than the others, with a flat base 
between ± 600 km s"^. There is also apparent emission at +900 km s~^ relative to systemic for 
Si XIV Lya (A6.182) . This is problematic since Si XIV Lya (A6.182) emission at +900 km s~^ 
would be difficult to reconcile with the small redshifts (consistent with systemic) of the [O Vll] 
and [Ne ix] emission features. All four profiles in Fig. 7 have widths that agree with the formal 
upperlimits in Table 1 (~ few hundred to ~ 1000 km s^^ FWHM for all profiles, except the more 
complex and broader Si xiv Lya (A6.182) profile which is ~ 2000 km s-i FWHM ). Formal results 
from Gaussian fits indicate that both O Vlll Lya (A18.969) and Mg Xll Lya (A8.421) are marginally 
resolved (see Table 1). 

Fig. 8 shows the velocity profiles of some of the strongest resonance absorption features and 
the profile of a transition in highly ionized Fe. Superimposed are velocity spectra from the best- 
fitting photoionization model (solid lines) and from the 99% upper limits on this model (dashed 
lines, superimposed in Fig. 8 (b,c) only) (sec §5 for further details). Fig. 8 (a) is centered on the 
Fe XX (A12.817) transition. Absorption is strong and spans ~ —300 km s^^ to ~ +500 km s~^. 
The profile appears deeper on the redward side of systemic. This may be due to a second absorber 
kinematic component, or due to blending with a one or more Fe L-shell transitions. Formally, from 
Gaussian fits, this Fe blend is marginally resolved (see Table 1). Fig. 8 (b) is centered on the 
Ne IX (r) Is^ — ls2p (A13.447) transition. The absorption trough spans ~ —100 km s~^ to ^ +300 
km s~^. Also noticeable in Fig. 8 (b) is the sharp absorption feature at ~ +1500 km s^^. This 
probably corresponds to a blend of Fe XVIII (A13.503) and Fe XIX (A13.504) L-shell transitions, 
which are at ~ +1300 km s^^ offset in velocity space from Ne IX (r) Is^ - ls2p (A13.4473). Fig. 8 

(c) shows the absorption velocity profile of Mg XI (r) Is^ — ls2p (A8.421). The profile minimum 
lies around ~ —230 km s~^, and the profile which spans ~ —600 km s"^ to systemic velocity, 
appears slightly broader on the blueward side, suggesting an unresolved kinematic component or 
an unresolved blend with (most likely) Fe XX L-shell (A9.166 — 9.167) transitions. Finally, Fig. 8 

(d) is centered on the Si xiii (r) Is^ — ls2p (A6.648) transition which spans ~ -300 km s-i to 
systemic velocity. 

Just as for the equivalent widths, one should not interpret apparent velocity shifts too literally 
because, in addition to uncertainty in the continuum, the velocity profiles are complex and in fact 
can be different for different ions. The complexity may either be intrinsic or due to contamination 
from other absorption and/or emission features. Emission filling in the red side of absorption 
features can lead to saturated features appearing non-saturated (Kaspi et al. 2002). If all the lines 
listed in Table 1 were in fact P Cygni, the estimated outflow velocities in Table 1 would gain a 
systematic error of about half of the listed outflow velocities. Intrinsic complexity could also be due, 
for example, to the presence of several kinematic components, making different contributions to 
the overall profile. However, the fact that most features are consistent with systemic (see Table 1) 
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is significant and those features that deviate strongly from systemic are Ukely affected by spectral 
complexity. 

5. Photoionization Modeling 

We used the photoionization code XSTAR 2.1.h to generate several grids of models of emission 
and absorption from photoionized gas in order to directly compare with the data. We used the 
default solar abundances in XSTAR (e.g. see Table 2 in Yaqoob et al. 2003). 

5.1. The Spectral Energy Distribution 

First we constructed a spectral energy distribution (SED) for NGC 4593 , as follows. Average 
radio, infrared (IR), optical and ultraviolet (UV) fluxes were obtained from Ward et al. (1987) and 
NED^^ . The X-ray region of the SED is constructed from the intrinsic continuum derived from our 
Chandra HETGS data (described below). Flux measurements obtained from historical data (Ward 
et al. 1987, NED) are not in general contemporaneous. Nevertheless, the NGC 4593 SED appears 
to exhibit a peak in the IR-UV spectral region. Due to the large spectral variability however, the 
flux levels in this spectral region can span up to an order of magnitude at a given frequency. Wc 
attempted to crudely account for variability in this region of the NGC 4593 SED by generating 
two different SEDs. For the first SED (SEDa), we chose the highest envelope of the measurements 
in the IR-UV region from Ward et al. (1987). SEDa exhibits a prominent IR-UV bump as is 
evident in Fig. 9 (SEDa is the solid line). We then simply connected the optical/UV endpoint of 
the IR-UV bump at ~ 10^'^ Hz to the X-ray portion of the SED with a straight line in log-log 
space. The second SED (SEDb) was generated simply by connecting the highest energy IR data 
point at ~ lO^^ *^ Hz from Ward et al. (1987) to the X-ray portion of the SED, such that SEDb does 
not have a prominent bump in the IR-UV region. SEDb is represented in Fig. 9 by filled circles 
connected by a dashed curve. By constructing these two SEDs, we could investigate the effect of 
reasonable uncertainties in the IR-UV region on the X-ray photoionization models. Also shown for 
comparison in Fig. 9 is the lower envelope of the measurements in the IR-UV region from Ward et 
al. (1987) (dash-dot line). 

To derive the X-ray portion of the SED, we fitted the data with a broken power law continuum 
and three absorption edges (as discussed in §3 above). The 0.5 keV intrinsic model flux was then 
simply joined onto the last UV (SEDa) and IR (SEDb) points of the lower-energy part of the SED 
by a straight line in log-log space. The hard X-ray power law was extended out to 500 keV. Note 
that we did not include any correction for the ACIS low energy CCD degradation as there is no 
satisfactory correction for grating data. Grids of photoionization models (see §5.2 below) were then 



^^NASA Extragalactic Database, http://nedwww.ipac.caltech.edu/ 
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generated using this SED. Also shown in Fig. 9 is the 'mean AGN' SED derived by Matthews k, 
Ferland (1987) normahzed to the ionizing luminosity of the NGC 4593 SEDa. 

5.2. Model Grids 

The photoionization model grids used here are two-dimensional, corresponding to a range in 
values of total Hydrogen column density, A'^h, and the ionization parameter, ^ = Lioa/ {rieB?)- Here 
Lion is the ionizing luminosity in the range 1 to 1000 Rydbergs, rig is the electron density and R 
is the distance of the illuminated gas from the ionizing source. The grids were computed for equi- 
spaced intervals in the logarithms of N^i and ^, in the ranges 10^0-22 ^m"^ and 10^-^~^ ° erg cm 
respectively. From SEDa, we deduce that Lion 7 x 10^^ ergs s~^ and therefore the lowest values 
of rig, log 1^ in the grid yield R ^ 150 pc and the highest values of rig, log in the grid yield R ^ 7 
light-hours. This range in R in the grid spans the expected location of the warm absorber. We 
computed grids with in the range 10^~^^ cm^^. For XSTAR models of the absorber, we confirmed 
that results from fitting the photoionization models to the X-ray data were indistinguishable for 
densities in the range rie = 10^ cm^^ to 10^^ cm~^. Hereafter we will use rie = 10^ cm~^ unless 
otherwise stated. Different density models will be important when we come to model the emission 
lines. 

All models assumed a velocity turbulence (6- value) of 1 km s~^, and since XSTAR broadens 
lines by the greater of the turbulent velocity and the thermal velocity, it is the latter which is 
relevant for these models. As the turbulent velocity increases, the equivalent width of an absorption 
line increases for a given column density of an ion. Therefore, to calculate model equivalent widths 
correctly one must know what the velocity width is (note that the 6- value could represent all sources 
of line-broadening). We cannot simply extend the XSTAR model grids into another dimension 
(velocity width), and deduce a 6- value directly from model-fitting, because we are limited by the 
finite internal energy resolution of XSTAR. The widths of the energy bins of the XSTAR model 
grids are in the range ~ 200-600 km s^^. Our choice of velocity width (6- value corresponding to 
the thermal width) ensures that the energy resolution of XSTAR is the limiting factor and not 
the 6-value itself. We will then need to take a rather more complicated approach to model-fitting 
(described in §5.3 below) and deduce a b-value appropriate for the data. 

5.3. Model-Fitting 

We proceeded to fit the MEG and HEG spectra simultaneously (binned at 0.04A, which is 
approximately the upper limit on the FWHM of most of the strong absorption features in Table 1) 
, with all corresponding model parameters for the two instruments tied together, except for the 
relative normalizations. The intrinsic continuum was a broken power law (the hard photon index 
again fixed at 1.794) modified by absorption from photoionized gas (derived from SEDa, which has 
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a prominent bump in the IR-UV region, see §5.1) and the Galaxy, as well as an absorption edge 
due to neutral Fe in the rest-frame of NGC 4593 . Excluding overall normalizations, there were a 
total of seven parameters, namely the column density of the warm absorber, Nh, the ionization 
parameter, ^ and the redshift z of the warm absorber, the energy E and maximum optical depth r 
of the neutral Fe absorption edge, the intrinsic soft photon index, r2 and the break energy of the 
broken power- law continuum, E-q- 

Our model-fitting strategy is more complicated than simply comparing the XSTAR model 
spectra with the data since the XSTAR spectra assume a certain line width (6- value), which itself 
is limited by the finite internal energy resolution of XSTAR. We chose 6=1 km in order that 
the energy resolution of XSTAR is the limiting factor. The result is that the equivalent widths 
of lines in the XSTAR spectra need to be calculated more rigorously using a line width (6-valuc) 
appropriate for the data. On the other hand, the ionic columns in the XSTAR models are robust 
from this point-of-view and do not depend on the velocity width. Therefore, after finding the best- 
fitting ionization parameter and warm-aborber column density from a global fit using the model 
described above, in § 5.5 we will deduce a 6-value from a curve-of-growth analysis, using ionic 
column densities from the best-fitting global model. We will then construct detailed model profiles 
for each of the detected absorption lines. 

The best-fitting parameters derived from the XSTAR global model fits were Nh = 5.37^1' ^1 x 
10^^ cm~^ and log^ = 2.52^q'q^ ergs cm s~^. The broken power-law model corresponding to 
this best fit had a break energy of E-q = I.OT^qq^ keV and a soft power-law photon index of 
Fi = 2.27^0 07, values which are consistent with the SED that was used as input to the models. 
The absorption edge had an energy of £^ = 0.708 ±0.002 and a optical depth of r = 0.19 ±0.04. The 
photoionized absorber was offset by —140 ± 35 km s~^ from the systemic velocity of NGC 4593 , 
the weighted mean of the values for velocity offsets of absorption features in Table 1, excluding the 
Fe XX (A12.817) transition, which is skewed, likely due to blending. We found that a photoionization 
model based on SEDb (the SED without a prominent bump in Fig. 9 ) yielded best-fit parameters 
within the 90% errors given above (AC ~ 4 for no additional parameters). Only a small difference 
in results between SEDa and SEDb should be expected since the relevant ionizing continuum is 
mainly in the X-ray band. Since there are no significant differences between the photoionization 
model fits corresponding to each SED, hereafter all photoionization model fits shall refer to SEDa 
(which has a prominent bump in the optical/UV frequency range). Table 2 shows the ionic column 
densities predicted by our best-fitting model for a range of ions. 

Fig. 10 shows the best-fitting model folded through the MEG response and overlaid onto the 
MEG counts spectrum. The lower panel of Fig. 10 shows the ratio of the data to the best-fitting 
model. Fig. 11 shows the actual best-fitting XSTAR model before folding through the instrument 
response. Fig. 12 shows the 68%, 90%, and 99% joint confidence contours of log^ versus Nh- 
Overall, the fit is good and it is already apparent that the best-fitting parameters which describe 
the overall spectrum also give good fits to the most prominent absorption lines and edges. This 
can also be seen in Fig. 3 which shows the best-fitting model overlaid on the photon spectrum, and 
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in Fig. 7 and Fig. 8 which show close-ups of the data and model (in velocity space), centered on 
particular atomic transitions. In Fig. 3 , some absorption features are discrepant with the model 
predictions, most notably Si Xlll (r) Is^ — ls2p (A6.648). The data and model disagree here for two 
reasons; firstly, because the 5-value used by XSTAR is not appropriate for the individual absorption 
features, as discussed above. Secondly, the data in Fig. 3 are binned at O.OSA which has combined 
the two minima in Fig. 7 (d) at ~ —200 km s"^ (which we identify with Si Xlll (r) Is^ — ls2p 
(A6.648) and ~ —1000 km (unidentified) into a single deep feature. 

The wavelength accuracy of the XSTAR atomic database is 5 mA, which corresponds to a 
velocity offset of 120 km s~^ at 1 keV. This offset is marginally less than the offset due to uncertainty 
in the MEG energy scale (5.5mA). We mentioned in §3 that the 2.0-2.5 keV region in the Chandra 
spectra suffers from systematics as large as ~ 20% in the effective area due to limitations in the 
calibration of the X-ray telescope absorption edges in that region We therefore omitted the 
2.0-2.5 keV regions of the Chandra spectra during the spectral fitting so that the fits would not be 
unduly biased. These and other instrumental features can be seen in the upper panel of Fig. 10 
. All the MEG data in the 0.5 5 kcV band are shown here, although the 2.0-2.5 keV region was 
omitted during spectral fitting. Other regions which have sharp changes in the effective area also 
show notable mismatches between data and model, but these are not as bad as in the 2.0-2.5 keV 
region. 

The ACTS CCDs have been undergoing a continuous low-energy QE degradation, probably 
due to absorption by a cumulative build-up of contaminants. Unfortunately, at the time of writing, 
there is no correction available for grating data. There is however, a model of the contamination 
pertinent for pure ACIS data (without the gratings). We can use this model as a 'worst case' 
to estimate the effects on our photoionization modeling and conclusions. Thus, we repeated the 
spectral fitting described above with the addition of the ACIS QE degradation model in XSPEC 
vll.2 (acisabs), using the defaults for absorption by the contaminants C, H, O, and N in the 
ratio 10:20:2:1 by number of atoms. Remarkably, we found that the best-fitting Nu, log^ and 
their associated statistical errors, were essentially indistinguishable from those obtained without 
the contamination model. What changed was the soft photon index, which increased to 3.04j;Q Q5. 
Even the break energy remained essentially the same. Thus, our photoionization modeling is driven 
by the discrete, narrow absorption features in NGC 4593 and is therefore robust to broadband 
calibration uncertainties. To determine why the degradation model did not change Nh or ^, we 
simulated data consisting of an absorbed power-law continuum with and without degradation 
respectively. We then fitted the data with and without the acisabs model. We found that when 
the simulated spectrum with no degradation was fitted with acisabs, Nh did not increase, but the 
continuum steepened. This is because the degradation yields deep absorption edges which are not 
in the simulated data. We conclude that the real data have much less degradation than is assumed 
by the acisabs degradation model. 



'http://asc.harvard.edu/udocs/docs/POG/MPOG/nodel3.html 
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5.4. Curve-of-Growth Analysis 

In order to determine the correct velocity width (6- value) of absorption lines for our models, we 
calculated curves-of-growth for different values of b (see Fig. 13 ). The theoretical curves-of-growth 
in Fig. 13 indicate that the measured equivalent widths of ~ lO^^'^A are most consistent with 
b ~ 100 — 200 km s~^. We also plotted the measured equivalent widths from Table 1 (converted to 
Angstroms) as a function of ionic column density (from the best-fitting XSTAR model described 
above) and overlaid these on the curvcs-of-growth (Fig. 13 ). For solar abundances, most of the 
measurements are consistent with b ~ 200 km s^^ within errors. Note that the predicted equivalent 
widths for b = 200 km s~^ are shown in Table 1 so that they can be compared directly with the 
observed equivalent widths. 

All of the plotted features in Fig. 13 are consistent with the curves of growth except for O Vll 
(r) Is^ - ls2p (A21.602). The disagreement with the curves of growth may be indicating a second 
component for the warm absorber with lower ionization parameter. However, all of the other strong 
absorption features on Fig. 13 seem to be consistent with the curves of growth. One explanation 
for the discrepancy might be an absorber component consisting of highly ionized dust. From §3 
we have seen that the data are consistent with an overabundance of neutral Fe, likely bound up 
in dust along the line-of-sight. The dust that lies closest to the AGN may well be highly ionized, 
contributing to a larger O Vll column than expected from a simple single zone photoionization 
model. Certainly this could also explain several highly ionized Fc features in Fig. 3 that arc not 
well fit by our best-fitting model in §5.3. Highly ionized dust can survive in warm, photoionized gas 
(see Reynolds et al. (1997) for further dicussion) although it seems probable in such scenarios that 
the warm material originates in the dusty cold material, which could be derived from a molecular 
torus. 



5.5. Detailed Comparison of Data and Photoionization Model 

In this section we compare in detail the best-fitting XSTAR photoionized absorber model and 
the data (refer to Fig. 3 for the spectrum, and to Fig. 7 and Fig. 8 for some velocity spectra). Note 

that although the best-fitting model was derived using spectra binned at 0.04A, we can compare 
the model and data at higher spectral resolution, depending on the signal-to-noise of the feature in 
question. The strongest absorption- line features in the Chandra data arc listed in Table 1. Fig. 3 
shows that the best-fitting XSTAR model described above {Nh = 5.37 x 10^^ cm^^, log^ = 2.52 ergs 
cm s~^ ) appears to underpredict the strength of several absorption lines. However, as outlined 
in §5.3, the equivalent widths of lines in the XSTAR spectra do not correspond to a (6-value) 
appropriate for the data. For example, Si xiii (r) Is^ — ls2p (A6.648) in Fig. 3 shows a large 
discrepancy between the model and the data, but Table 1 does not. A proper comparison between 
data and model for the stronger lines is described below. 

The model velocity profiles in Fig. 7 (for four of the strongest absorption lines: O Vlll Lya 
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(A18.969) , Ne x Lya (A12.134) , Mg Xll Lya (A8.421) & Si XIV Lya (A6.182) were calculated using 
the appropriate 6-valuc for the data (200 km s~^) and show reasonable agreement between data 
and model. In this case the model curves were calculated using the broken power-law and two-edge 
continuum described in §3, a Gaussian for the absorption line, with a = b/V2 = 141.4 km s~^, 
and the predicted model equivalent width from Table 1. Fig. 8 , which shows the velocity profiles 
of some of the strongest resonance absorption features and the profile of a transition in highly 
ionized Fe, illustrates that the best-fitting XSTAR model under-predicts some features. However 
the lower limits on the best-fitting column density and ionization parameter at 99% confidence, 
yield ionic columns that give velocity profiles that are more consistent with the data for Nc IX 
(r) Is^ - ls2p (A13.447) and Mg xi (r) is^ - ls2p (A8.421) (dashed lines in Fig. 8 ), without 
significant disimprovement of the profiles in Fig. 7 . Thus, in spite of the likely wide range of 
ionization states in the warm absorbing gas, a single-component photoionization model can give 
surprisingly good fits to the principal absorption lines. 

There are many more discrete absorption features apparent in the spectra but with lower 

signal-to-noise. Although these individual weak features cannot be studied in detail we can at least 
check whether the best-fitting XSTAR model does not conflict with the data. We found that, as 
far as narrow absorption lines are concerned, most absorption lines predicted by the model were 
consistent with the data. In particular, we found that the model agrees well with the data for the 
low-order Lyman and He- like resonance transitions of Ne and O, although there may be blending 
associated with higher order transitions such as Ne IX Is^ — ls4p (All. 000). Conversely, the data 
show possible strong absorption features which are not predicted by the model. These features are 
mainly in the ~ 1.1 — 1.5 keV range (e.g. see Fig. 3 (d)). Some of the strongest features may 
correspond to blends of L-shell transitions in very highly ionized Fe (Fe XX -Fe XXV ). The atomic 
database used by XSTAR 2.1.h does not yet include all transitions in highly ionized Fe and this 
would certainly explain these discrepancies between data and model. 

The Fe ionic columns predicted by the best-fitting XSTAR model are largest for ionization 
states Fe XIX to Fe XXIV . According to Behar, Sako, & Khan (2001), an Fe M-shell unresolved 
transition array (UTA) due to Fe l-Fe XV lies between 17.5A and 15A. According to Gu et al. 
(2001), there is a strong L-shell signature due to Fe xx - Fe XXIV between 10. SA and 12. 5A. 
Thus, some of the strong absorption features between 0.85-1.0 keV in Fig. 3 may be part of, or 
blended with, an Fe UTA due to highly ionized Fe (our Fe ionic columns are largest for Fe XIX - 
Fe XXV ). UTAs due to less highly ionized Fe (in the ~ 15-17 A range) have been detected in some 
other Seyfert 1 galaxies (IRAS 13349-^2438, Sako et al. 2001a; MCG -6-30-15, Lee et al. 2001; 
NGC 3783, Kaspi et al. 2002; NGG 5548, Kaastra et al. 2002; Mkn 509, Yaqoob et al. 2003). 

In §3 we mentioned how spectral fits using simple absorption-edge models (representing O Vll 
and O Vlll) could be misleading. We see that more sophisticated modeling with better data also 
gives a weak O vii edge, as obtained from edge-fitting, but also a much weaker O viii edge since 
the data above the O Vlll edge are affected by other absorption features (Fe UTAs and Ne IX (r) 
(A13.447) in particular) which can be accounted for by the photoionization models. This has also 
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been pointed out by Lee et al. (2001). The complex of O Vll Is^ — Isnp resonance transitions 
with n > 5 may also complicate the region just above the O Vll edge (see also Lee et al. 2001). 
From the XSTAR models wc can extract the O vii and O viii ionic column densities at the 
extremes of the 90% confidence intervals for Nh versus ^. We find that the O vii column is 
4.01^ 5 ^ 10^^ cm~^ and the O Vlll column is 2.7^^^ x 10^^ cm~^. Given an absorption cross- 
section at threshold of 2.4 x 10~^^ cm~^ for O Vll and 9.9 x 10~^° cm~^ for O Vlll (Verner et al. , 
1996), we find that the model optical depths (at threshold) of the O edges are tq vii — '-'■001^0 ooo5 
and To viii = 0.026lQ'Qgy respectively. These values can be compared with those obtained from an 
ASCA observation of NGC 4593 in 1994, in which simple 'edge-models' yielded tq vii = 0.26 ± 0.04 
and TO viii = 0.09t|j;^^ ((Reynolds 1997); see also George et al. 1998). Whilst the ASCA result is 
inconsistent with our inferred O Vll optical depth, it is consistent with the optical depth of our best 
fit neutral Fe L edge (r = 0.19 it 0.04). Note that ASCA resolution at energies around 0.7 keV is 
too coarse to resolve the two edges separately (see discussion in §3). From an Fe L3 edge depth of 
0.19 lb 0.04, and given an Fe 2p absorption cross-section of 1.25 x lO^^^cm^^ (Verner et al. 1996), 
wc infer a neutral Fe column of 1.5ib0.3 x 10^'' cm^^. An Fe LI (2s) edge (rest energy of 0.845 keV) 
optical depth of r < 0.05 is allowed by the data, implying a neutral Fe column of < 3.6 x lO^^cm^^ 
which is consistent with the estimate from the Fe L3 edge, although we caution that the spectrum 
at both the Fe LI and the Fe L2 (rest energy of 0.720 keV) edge energies is complex. Indeed, we 
cannot derive a reasonable constraint on the Fe L2 edge optical depth because of apparent emission 
blueward of this edge (see Fig. 3 ). The implied neutral column density for cosmic abundances is 
then Nh ~ 3.8 x 10^^ cm~^, which would act to extinguish the soft X-ray spectrum. Therefore the 
neutral Fe could take the form of neutral dust, as argued by Lee et al. (2001). Our estimate of the 
O VIII optical depth from the photoionization modelling is inconsistent with both the ASCA result 
and our direct fit of the O Vlll edge (r = 0.14+o;q3). However, the XSTAR model fits the data well, 
so simple edge models may give the wrong optical depth. Nevertheless, we find that the measured 
equivalent widths of the O Vll (r) Is^ - ls2p (A21.602) and O Vlll Lya (A18.969) absorption fines 
in Table 1 yield ionic column densities that are consistent with the edge optical depths measured 
in §3. 

We investigated the robustness of the XSTAR model fits to the absorption features by examin- 
ing the data versus model at the extremes of the Nh versus ( confidence intervals (see Fig. 12 ). By 
picking pairs of values of Nh and ^ in the plane of the contours, we can freeze the column density 
and ionization parameter at the chosen values and find the new best-fit for each pair of values (this 
is in fact how such contours are constructed). Then we can examine the data and model spectra 
in detail for each of these new fits to see where the model and data deviate significantly from each 
other. We found that the models were still able to account for the strengths of all of the strong 
absorption-line features (listed in Table 1) for pairs of Nh and ^ lying inside the 99% confidence 
contour. On the other hand, for pairs of Nh and ^ lying outside the 99% contour we found that 
absorption lines in the 12-15A range were particularly sensitive and the equivalent width would be 
over-predicted or under-predicted by as much as a factor of two. Therefore, we conclude that the 
99% confidence contours in Fig. 12 are a good guide to the allowed range of Nh and ^. 
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As well as absorption, the XSTAR code also calculates emission spectra. However, the very 
weak [Nc ix] and [O Vll] forbidden lines that wc observe (see Fig. 6 ) cannot constrain the XSTAR 
emission spectra. Therefore our best-fit photoionization model of the warm absorbing gas (Fig. 11 
) does not warrant inclusion of the corresponding (weak) emission spectrum. Stronger emission 
features have been detected in other observations of NGC 4593 and we shall discuss potential 
constraints on the emission region from previous observations in §6 below. The possibility of very 
broad relativistic emission lines in our soft X-ray data shall also be discussed in detail in §7. 

6. Comparison with previous observations 

Over the past decade NGC 4593 has been observed in X-rays, by ASCA on 1994 January 9, 
BeppoSAX from 1997 December 31 to 1998 January 2 and most recently, by XMM-Newton and 
Chandra LETGS (non-simultaneously on 2000 July 2 and 2001 February 16 respectively). We 
measured a 2-10 keV flux of 4.0 x 10"^^ ergs cm^^ g-^ from NGC 4593 with Chandra HETGS 
. This is comparable with the flux detected by BeppoSAX (Guainazzi et al. 1999) and ASCA 
(Reynolds 1997). Using Fig.l of Kaastra k. Steenbrugge (2001), we find compared to the Chandra 
LETCS and XMM-Newton observations, our MEG flux at 20l is a factor ~ 1.6 and ~ 1.9 higher 
respectively. We compared our soft X-ray measurements with those from the ASCA observation 
in §3 and shall not discuss the ASCA results further. BeppoSAX observed a similar underlying 
continuum to that observed by us with Chandra HETCS . When the BeppoSAX continuum is fit 
with a simple powerlaw model (F = 1.70 ±0.04), the main deviations include an absorption feature 
at ~ 0.7 keV and an apparent soft excess between ~ 0.3-0.7 keV. Guainazzi et al. (1999) measured 
two edges whose threshold energies and optical depths are consistent with the Chandra Fe L and 
O VII measurements. 

A strong soft excess was also observed by Chandra LETCS at energies less than ~ 0.65 keV, 
when the power law describing the LETCS data at 1 — lOA is extrapolated to lower energies. 
The Chandra LETGS data also exhibit a flux deficit between ~ 0.65 — 1.2 keV relative to the 
extrapolated powerlaw. Kaastra & Steenbrugge (2001) claim that the Chandra LETCS spectrum 
exhibits few strong discrete emission or absorption features, unlike our data. We note from Fig. 1 
in Kaastra k, Steenbrugge (2001) (which shows counts spectra versus observed wavelength between 
18-23 A), that the Chandra LETCS spectrum docs exhibit some absorption features. There is 
weak absorption due to O vii (r) Is^ — ls2p (A21.602), a weak [O vii] emission feature and some 
absorption due to highly ionized Fe (possibly a blend of Fe xxiii (A21.400) and Fe xxiv (A21.383) 
inner shell transitions). Lower statistics during the Chandra LETCS observation may explain the 
weakness of many such features. In particular, Kaastra &; Steenbrugge (2001) note the weakness of 
any possible absorption lines from highly ionized C, N and O. They find from photoionized models 
that ionic columns of the H- and He-like ions of C, N, O are < 10^® cm~^. These ionic column 
densities are significantly lower than the 10^^ — 2 x lO^'' cm~^ range we find for H-likc C, N, O ions 
from our best-fit model in §5.3. Kaastra &: Steenbrugge (2001) report that the best-fit model for 
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the warm absorber predicts continuum absorption edges of less than a few percent for O Vll and 
O VIII . This agrees with our upper hmit on the O VII edge depth of < 0.03, but is inconsistent 
with our estimate of the O Vlll edge depth (0.14^q'q4). Kaastra & Steenbruggc (2001) argue that 
one explanation of the flux reduction in the 12-15A band in the Chandra LETGS data is a highly 
ionized warm absorber with Nh ~ 2.4 x 10^^ cm~^ and log^ ~ 2.54 ergs cm s~^. This value of 
the ionization parameter is in excellent agreement with that of our bestfit model, although their 
column density is less than half the value we observe. The lack of strong edges in the LETGS 
and XMM-Newton RGS data, and the much lower column density of their preliminary best-fitting 
warm absorber model, according to Kaastra &: Steenbruggc (2001), suggest that significant changes 
occurred in the warm absorber in the ~ 230 days between the LETGS and HETGS observations. 

The details of the continuum observed and spectral features detected during the short observa- 
tion of NGC 4593 with XMM-Newton RGS are not discussed by Kaastra k Steenbruggc (2001). We 
note however, from Fig.l in Kaastra & Steenbruggc (2001), that the RGSl data appear to exhibit 
a relatively strong [O Vll] emission feature, unlike the Chandra LETGS data, and there are hints 
of a P Cygni profile of O Vlli Lya (A18.969) . Also from Fig.l in Kaastra & Steenbruggc (2001), 
there is weak absorption in the RGSl data due to O vii (r) Is^ — ls2p (A21.602) and highly ionized 
Fe. The presence of the [O vii] emission line in the XMM-Newton RGSl data and its absence in 
the Chandra LETGS data from a year later, suggests that the O Vii emission region lies within 1 
light year of the source of the continuum radiation. The [O Vll] emission feature likely originates in 
a photoionized plasma and the reduction in its intensity over the year between the XMM-Newton 
RGS observation and the Chandra LETGS observation is consistent with the halving of the mean 
flux from this source in that time. 



7. Dust or Relativistically Broadened Lines? 

In the MEG spectrum of NGC 4593 we observe a sudden drop in continuum flux of 20% at 
~ 0.7 keV (apparent in Figs. 1, 2 and 3) which we identify with an Fe L3 edge. A similar edge-like 
feature has been observed at the same (rest-frame) energy in MCG —6—30—15 and Mrk 766 with 
XMM-Newton RGS (Branduardi-Raymont et al. 2001 (hereafter BROl), Sako et al. 2001b, Mason 
et al. 2002) and independently with Chandra HETGS (Lee et al. 2001). BROl, Sako et al. (2001b) 
and Mason et al. (2002) argue that the edge-like feature in particular, and more generally the soft 
excess in some Seyfert 1 galaxies, can be explained by relativistically-broadcncd emission lines due 
to C VI (A33.736), N vii Lya (A24.781) and O viii Lya (A18.969) . Lee et al. (2001) argue that 
the edge-like feature is due to a combination of neutral Fe bound-free L-shell absorption edge due 
to neutral dust and higher order (n > 5) O Vll resonance Is^ — Isnp absorption. In this section 
we shall consider our observations in the context of each interpretation. 
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7.1. Testing the relativistic emission model for NGC 4593 

Due to the MEG bandpass, we can only test for relativistically-broadened O viii Lya (A18.969) 
emission in our NGC 4593 data. Firstly, we added an emission-line model from a relativistic disk 
around a Schwarzschild black hole (e.g. Fabian et al. 1989) to the best-fitting photoionized absorber 
model (with a broken power law model of the continuum but vjith the Fe L3 edge removed). The 
inner and outer radii of the disk were fixed at 6 and 1000 gravitational radii respectively, and the 
radial line emissivity per unit area was a power law with index —2.5. The line energy was fixed at 
0.6536 keV in the NGC 4593 rest frame. The best-fitting inclination angle and equivalent width 
were 39 ± 1 degrees and ISlg eV respectively. The blue wing of the Schwarzschild disk line exhibits 
a sharp drop and gives a good overall fit to the sharp drop in flux at ~ 0.7 keV (see Fig. 3 ), 
although the apparent flux deficit above ~ 0.7 keV in Fig. 3 is somewhat underpredicted. The 
C-statistic of the best-fit Schwarzschild disk line was worse than that for the best-fit Fe L3 edge 
model by ~ 30 for the same number of parameters. We note that if a blackbody is used for the 
soft continuum instead of a powerlaw, we obtain the same inclination angle but a larger equivalent 
width (24^7 eV) and a C-statistic which is still larger by 29 than the Fe-L edge model. 

When we replaced the O viii Lya (A18.969) emission from around a Schwarzschild black 
hole with the corresponding emission from around a maximally spinning Kerr hole, we found that 
the C-statistic worsened by ~ 6 relative to the Schwarzschild disk line, for the same number of 
parameters. The best fitting inclination angle and equivalent width were 771^0 de grees and 11 it 7 
eV respectively. We assumed the same radial emissivity law but extended the inner radius down 
to 1.24 gravitational radii. The Kerr disk line cannot fit the data around 0.7 keV because the blue 
wing of the Kerr disk line, being strongly broadened, is less sharp and 'edge-like' than that of the 
Schwarzschild disk line. 

A problem with the disk line interpretation is that if the drop in flux occurs at the same energy 
in several sources (and there are now three examples of such), the implied disk inclination angle 
would have to be identical for all the sources within 1°. An obvious test of the neutral dust model 
would be to measure the O I column predicted by the inferred Fe I column if it were locked up in 
dust. Unfortunately, due to contamination, the soft X-ray quantum efficiency of Chandra ACIS 
has degraded, and produces significant excess absorption due to neutral O. Neutral dust along the 
line-of-sight might also reasonably be expected to contain silicates. We tested the data for a neutral 
Si I edge and found that, with the current calibration, r < 0.06 is allowed, although this is not 
significantly detected. This yields a Si I column of < 4.4 x 10-*^^ cm~^ which is consistent with the 
inferred Fe i column of ~ 1.5 x 10^'' cm^^. Alternative mechanisms for obtaining an overabundance 
of Fe I along the line-of-sight to the AGN could include shattered stellar cores from an environment 
rich in supernova remnants or a progenitor gamma-ray burst associated with the host galaxy (Piro 
et al. 1999). 
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8. Summary of Observational Results 

We observed NGC 4593 for ~ 79 ks with the Chandra HETGS on 2001 June 29, simultaneously 
with RXTE . The complex Fe-K line and Compton-reflection continuum will be discussed elsewhere. 
Here we summarize the main results from the Chandra HETGS soft X-ray spectrum, in view of 
relevant constraints from RXTE . 



1. The 0.5-5.0 keV data can be modeled with a broken power-law intrinsic continuum plus a 
single zone warm absorber model (see below). RXTE was used to constrain the hard power- 
law index (1.794), yielding a soft X-ray index of 2.27lQ;Qy and a break energy of 1.071q;q7 
keV when self-consistently modeled with the photoionized absorber. During the observing 
campaign the 2-10 keV flux was 4.0 x 10~^^ ergs cm~^ s~^, typical of historical values (the 
corresponding 2-10 keV luminosity was 6.1 X 10^2 gi-gs g-i^ foj. ^ 70 km s'^ Mpc"^ and 
qo = 0). 

2. Below ~ 2 keV, the observed spectrum shows considerable deviations from a smooth con- 
tinuum, due mainly to bound-free absorption opacity, and complexes of discrete absorption 
features. Particularly noteworthy are the Fe L3 and O viii edges at 0.708 =b 0.003 keV and 
0.868to o?o respectively. Within the errors, both values are consistent with their respec- 
tive rest energies, and from the best-fitting edge models the optical depths at threshold are 
0.18 ± 0.04 for Fe L3 and 0.14l°;°| for O Vlll. The optical depth of the Fe L3 edge is consis- 
tent with the depth ascribed by ASCA to an O vil edge, whilst our O viii edge depth is in 
marginal agreement with the ASCA measurement ((Reynolds 1997); (George et al. 1998)). 
We note that ASCA spectral resolution was incapable of distinguishing an Fe L3 and O vii 
edge and that estimates of the O Vlll edge depth are also complicated by nearby absorption 
features due, for example, to Ne and Fe. The neutral Fe edge may be due to dust in the warm 
absorber and therefore we tested the data for neutral Si. We found that, with the current 
calibration, a neutral Si i edge with r < 0.06 is allowed, which yields a Si i column that 
is consistent with the inferred Fc I column. The soft X-ray quantum efficiency of Chandra 
ACIS has degraded due to contamination, and produces significant excess absorption due to 
neutral O, so we could not deduce whether there was excess absorption due to neutral O 
in NGC 4593 . An alternative model of the neutral Fe edge region due to relativistically 
broadened emission lines from an accretion disk yielded a significantly worse fit to the data 
than an Fe L edge. 

3. On the next level of detail, we detect strong discrete absorption lines from H- and He-like 
ions of O, Nc, Mg, Si as well as discrete absorption features due to highly ionized Fe XIX- 
XXV. Only the N vii Lya (A24.781) , O viii Lya (A18.969) and Mg xii Lya (A8.421) lines 
appear to be marginally resolved (see Table 1). A marginally resolved absorption feature at 
~ 12.82A is likely to be a blend of Fe XX lines. The remaining strong, discrete absorption 
features are unresolved (FWHM< 275 km s~^ for O Vll (r) Is^ — ls2p (A21.602), going up 
to FWHM< 2145 km s~^ for the complex Si XIV Lya (A6.182) absorption feature). The 
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absorption lines are consistent with a net outflow velocity of ~ — 140 km s^ relative to 

systemic. However, the detailed velocity profiles of the absorption features are not all the same 
for all ionic species. The complexity and differences in profiles could be due to the contribution 
from different, unresolved kinematic components and/or blending and contamination from 
absorption/emission due to other atomic transitions. 

4. Wc also detect O vii (r) Is^ - ls2p (A21.602) and N vii Lya (A24.781) in our local frame, 
possibly indicating absorption in our Galaxy or local group of galaxies due to a hot medium. 
This absorption is likely to be due to the same medium as that responsible for the local 
{z = 0) absorption observed by Nicastro et al. (2002). 

5. We have modeled the Chandra HETGS spectra using the photoionization code XSTAR 2.1.h. 
The best-fitting ionization parameter and equivalent Hydrogen column density of the absorb- 
ing gas are log^ = 2.521q;o4 ergs cm s~^ and Nh = S.STtJl^g x lO^^cm"^ respectively. This 
best-fitting model gives a good fit to the overall observed spectrum. It also is able to model 
most of the principal local absorption features well but underpredicts some absorption lines 
such as Ne IX (r) Is^ — ls2p (A13.447) by up to a factor of ~ 2 in equivalent width, but 
these are likely cases of blending with lines of highly ionized Fe. A curve-of- growth analysis 
indicates that a velocity width of 6 ~ 200 km s~^ is consistent with the data. Our results 
are insensitive to reasonable variations in the unobserved UV part of the spectral energy 
distribution. 

6. We detect weak emission features due to [O Vll] and [Ne ix] . However, our data cannot 
constrain the distance of the emitter from the ionizing source. Since we cannot assume that 
the emitter and absorber have the same ionization state and/or column density, time-resolved 
spectroscopy with better signal-to-noise is required to address the question of the location of 
the absorber and emitter. Deducing the mass outflow of the X-ray absorber must also await 
better data because the density and global covering factor are unknown. However, we obtain 
upper limits on the electron density of the emitting material of < 8 x 10^° cm~^ from the 
ratios of the [O Vll] and O Vll (i) line strengths. We note also from observations of NGC 4593 
with XMM-Newton and Chandra LETGS (see §6), that the [O Vll] line emission is variable 
on a timescale less than a year. Therefore the [O Vll] emission region appears to lie within 1 
light-year of the ionizing radiation source. 



9. Conclusions 

We observed the Seyfert 1 galaxy NGC 4593 with the Chandra high energy transmission 
gratings and measure strong absorption lines from He-like O, Ne, Mg, Si, H-like N, O, Ne, Mg, Si 
and highly ionized Fe XIX- XXV. The offset velocity of the strongest absorption profiles are slightly 
blueshifted but the errors are consistent with the systemic velocity of NGC 4593 . Many profiles 
hint at the presence of either multiple kinematic components or blending. Absorption features 
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due to N VII Lya (A24.781) , O Vlll Lya (A18.969) and Mg Xll Lya (A8.421) as well as a blend 

of Fe XX features appear to be marginally resolved. The soft X-ray spectrum of NGC 4593 is 
adequately described by a simple, single-zone photoionized absorber with an equivalent Hydrogen 
column density of 5.37^q^q x 10^^ cm~^ and an ionization parameter of log^ = 2.52^qq^ ergs 
cm s"^ although there remain some features which are not identified. Although the photoionized 
gas almost certainly is comprised of matter in more than one ionization state and may consist of 
several kinematic components, data with better signal-to-noise ratio and better spectral resolution 
are required to justify a more complex model. In emission we detect only weak forbidden lines of 
[Ne ix] and [O Vll] . We identify a spectral feature at ~ 0.707 keV with a neutral Fe L edge, which 
is most likely due to neutral dust along the line-of-sight to NGC 4593 . However, this is not the 
only interpretation of this feature and a search for neutral O absorption which would reasonably 
be expected from dust absorption is complicated by contamination of the Chandra acis CCDs. 
Neutral Si absorption is allowed by the data and corresponds to a neutral Si column density that 
is consistent with the inferred neutral Fe column. We detect N Vll Lya (A24.781) and O Vll (r) 
Is^ — ls2p (A21.602) absorption at z = 0, due to a local, hot medium. The hot, local [z = 0) 
medium responsible is likely the same hot medium responsible for the absorption seen by Nicastro 
et al. (2002). 
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Tablc 1. Spectral Lines in the Chandra HETGS Spectrimr of NGC 4593 
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Note. — Absorption-line parameters measured from the MEG spectrum, using simple Gaussians (see 
§4). All measured quantities refer to intrinsic parameters, already corrected for the instrument response. 
Errors are 90% confidence for one interesting parameter (AC = 2.706). All velocities have been rounded 
to the nearest 5 km s~^. " Laboratory- frame wavelengths. Measured equivalent widths in the NGC 4593 
frame. ^ Predicted equivalent widths using b = 200 km s-i and XSTAR columns (§5.5). Velocity offset 
(NGC 4593 frame) of Gaussian centroid relative to systemic. Negative values are blueshifts. Change in 
the fit statistic on addition of Gaussian. ^ Velocity offset (galactic frame) of Gaussian centroid relative to 
systemic. ^ Emission line. 
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Table 2. Ionic Column Densities (10^^ cm ^) 
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IQi'^ predicted by the best-fitting XSTAR 

model to the Chandra data (see §5). The pre- 
dicted column densities for Si iii-iv and C ii-iii 
from the XSTAR model were negligible. 
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Figure Captions 

Figure 1 

The NGC 4593 MEG spectrum binned at 0.32A compared with the best-fitting power-law model 
from a joint fit of HETGS and 3-8 keV RXTE data, extrapolated down to 0.5 keV. 

Figure 2 

Comparison of the NGC 4593 MEG data (binned at 0.32A) to a broken power-law model modified 
by two absorption edges and Galactic absorption (1.97x 10^° cm~^). 'Edge-models' have historically 

been used to model lower spectral resolution CCD data. The hard X-ray photon index is fixed at 
1.794, and the best-fitting soft X-ray index and break energy are 2.26 it 0.05 and LOS^q q^ keV 
respectively. The best-fitting edge energies of 0.708± 0.003 keV and O.SGSIqqJq keV agree well 
with the expected rest-frame energies of the neutral Fe L3 (0.707 keV) and O Vlll (0.871 keV) 
edges respectively. The addition of an O Vll edge (0.739 keV) does not significantly improve the 
fit. The threshold optical depth of the O Vll edge is < 0.03 at the 90% confidence level. Note 
that using this model, the inferred optical depth of the O Vlll edge could be affected by the Nc IX 
Is^ — ls2p (A13.447) absorption feature (at ~ 0.89 kcV, observed). The inferred depths of the 
O VII and Fe L3 edges could also be influenced by complexity in the spectrum between ~ 0.7 — 0.8 
keV, in particular from Fe inner-shell and O Vll resonance absorption. 

Figure 3 

NGC 4593 MEG observed photon spectrum (binned at O.OSA) compared to the the best-fitting 
photoionized absorber model (red solid line). Also shown is the intrinsic continuum (blue solid 
line) modified by Galactic absorption (neither the data nor model have been corrected for Galactic 

absorption). The dashed lines show the expected positions of some bound-free absorption edges. 
The model consists of an intrinsic continuum which is a broken power law (best-fitting break energy 
at 1.07 keV), absorbed by photoionized gas with best-fitting ionization parameter of log^ = 2.52 
ergs cm and column density 5.37 x 10^^ cm~^. Full details of the model calculations, fitting 
procedures, and discussion of the details of the comparison between data and model can be found 
in §5. The 0.7-1.1 keV region is very complex so simple two-edge models fitted to older CCD 
data (of Seyfert 1 galaxies in general) could have been biased by this complexity. There is also an 
instrumental feature between ~ 0.85-0.9 keV (see Fig. 10 ) which may partly be responsible for 
the poor fit in this region. 

Figure 4 

Chandra MEG photon spectrum for NGC 4593 against wavelength in the source rest-frame, using 
z = 0.0083. The data are binned at 0.02A, approximately the MEG FWHM spectral resolution 
(O.O23A). Labels show the Lyman series wavelengths (in blue) for Ar, S, Si, Mg, Na, Ne, O, and 
N. Also shown (in red) are the wavelengths of the Helium-like triplets (resonance, intercombination 
and forbidden lines) of Ar, S, Si, Mg, Ne and O. 

Figure 5 

Same spectrum as Fig. 4, but shown in blue are the wavelengths of He-like resonance-absorption 
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transitions (n = 1 ^ 2,3,4,..) of Ar, S, Si, Mg, Na, Ne, O and N. Also shown (in red) are K-shell 
transitions of Fe XIX in the wavelength range 2 — 25A, as an illustration of some of the Fe blending 
in the data. 

Figure 6 

Chandra MEG photon spectrum for NGC 4593 against observed wavelength in the vicinity of 
the He-like Ne and O triplets (binned at 0.04A and 0.02A respectively). Plotted in red are the 
wavelengths of the respective He-like resonance, intercombination and forbidden transitions. 

Figure 7 

Velocity profiles from combined Chandra MEG and HEG spectra (except for O viii Lya (A18.969) 
which is from MEG data only since it is outside of the bandpass of the HEG). The profiles correspond 
to some of the strongest absorption features in the data (see Table 1). The absorption features are 
well-described by the column densities from the best-fitting XSTAR photoionization model (black 
solid lines), which has log^ = 2.52 ergs cm s~^ and Nh = 5.37 x 10^"*^ cm~^ (see §5). The model 
profiles were calculated using the model EW and a Gaussian with a = 6/ a/2, for b = 200 km s~^, 
inferred from a curve-of-growth analysis (see §5.5 and Fig. 13 ). The profiles were given a uniform 
blueshift of —140 km s~^, corresponding to the weighted mean of the velocity shifts in Table 1, 
excluding the Fe XX absorption feature which is skewed due to likely blending (see §5.5). A velocity 
of zero corresponds to the systemic velocity for NGC 4593 (assuming z = 0.0083). A negative 
velocity here indicates a blueshift relative to systemic. The centroids of the absorption features 
lie between systemic and ~ —350 km s~^ (see Table 1), but all the profiles extend redward to the 
systemic velocity and slightly beyond. FWHM velocity widths from Gaussian fitting are given in 
Table 1 and can be compared with the FWHM MEG velocity resolution of 360, 565, 810 and 1105 
km s-i for O Vlll Lya (A18.969) , Ne X Lya (A12.134) , Mg Xll Lya (A8.421) and Si XIV Lya 
(A6.182) respectively, at observed wavelengths. 

Figure 8 

As Fig.7, for four more velocity spectra from combined Chandra MEG and HEG data. Solid 
lines correspond to absorption profiles calculated from the best-fitting photoionization model as in 
Fig. 7 . Dashed lines correspond to absorption profiles calculated using the upper limits (at 99% 
confidence) of the column density from the best-fitting photoionization model. Clearly the best-fit 
photoionization model under-predicts some of the features which the spectra are centered on (as 
labeled) but the upper limits on the ionic column densities are more consistent with the profiles for 
Ne IX (r) Is^ - ls2p (A13.447) and Mg xi (r) Is^ - ls2p (A9.169) in particular. The profiles were 
given a uniform blueshift of —140 km s^^, which corresponds to the weighted mean of the velocity 
offset of the strongest absorption features (see §5.5). FWHM velocities from Gaussian fitting are 
given in Table 1 and can be compared with the FWHM MEG velocity resolution of 535, 510, 745 
and 1030 km s-^ for Fe XX (A12.817), Ne IX Is^ - ls2p (A13.447), Mg XI Is^ - ls2p (A9.169) and 
Si XIII Is^ — ls2p (A6.648) respectively, at observed wavelengths. Blending with Fe transitions 
is particularly likely to be associated with the Fe XX (A12.817) and Ne IX Is^ - ls2p (A13.447) 
transitions. 
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Figure 9 

The observed and constructed baseline SED (solid curve) used for photoionization modeling of 
NGC 4593 . Average radio and IR fluxes arc from Ward et al. (1987). The optical and UV data 
are fluxes from non-contemporaneous observations in NED. The X-ray portion is constructed from 
our Chandra HETGS and RXTE data. The hard X-ray power law has V = 1.794 and has been 
extended to 500 keV. Given the significant variability in the UV/optical we have constructed upper 
and lower envelopes enclosing the data. The solid and dash-dot lines correspond to the upper 
(SEDa) and lower envelopes of optical-UV data in NED respectively. Note the prominent bump in 
the IR-UV region of the spectrum. The dashed line (SEDb) corresponds to simply joining the IR 
flux at 0.33 eV and the X-ray flux at 0.5 keV with a straight line in log-log space and thus has no 
prominent bump in the IR-UV region. See §5.1 for full details of the construction and applications 
of the SEDs. The dotted curve is the 'mean AGN' SED of Matthews & Ferland (1987), normalized 
to the same ionizing luminosity (i.e. in the range 1-1000 Rydberg) as that of the NGC 4593 baseline 
SED (SEDa). 

Figure 10 

Comparison of MEG data (top panel, raw counts spectrum) binned at 0.16A with the best-fitting 
XSTAR photoionization model (red curve; see §5) and the ratio of data to model (bottom panel), 
showing a good overall fit to the data. In this representation the various instrumental edges and 
features can be seen, and in particular the poor fit around 2.0-2.5 keV can be seen to coincide with 
the large jump in the X-ray telescope effective area. 

Figure 11 

Best-fitting XSTAR (v.2.1.h) photoionization model spectrum fitted to the MEG and HEG data 
as described in §5 (see also Fig. 3 ). The ionization parameter is log^ = 2.52 ergs cm s~^, and 
Nh = 5.37 X 10^^ cm~^. The intrinsic continuum is a broken power law. A fitted neutral Fe edge at 
0.708 keV, and Galactic absorption are also included. The purpose of this plot is to show the effects 
of the instrumental effective area, especially where it changes abruptly. Note that this spectrum is 
for b =1 km s~^ so the equivalent widths of absorption features are limited by the energy resolution 
of XSTAR and are not as accurate as those computed for Fig. 7 and Fig. 8 (to which the reader 
should refer for details of the absorption features). Full details and other parameters can be found 
in the text (in §3 and §5). 

Figure 12 

The 68%, 90%, and 99% confidence contours of the logarithm of the ionization parameter (log^) 
versus the column density (Nh) around the best-fitting photoionized absorber corresponding to 
Fig. 10 (see also Fig. 3 and Fig. 11 ). The HEG and MEG spectra were binned at 0.041. See §5 
for details. 

Figure 13 

Curves of growth for various values of the velocity width, b. Plotted for each b is the logarithm 
of line equivalent width (EW) per unit wavelength against the logarithm of the product of ionic 



-31 - 



column density (N), oscillator strength (/), and wavelength (A). The plotted points are taken from 
the measured equivalent widths (converted to A using the values in Table 1) and predicted ionic 
column densities (from the best-fitting XSTAR model described in §5). These values are consistent 
with b =200 km within errors. 
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